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Two-sided flux decoration experiments indicate that threading dislocation lines (TDLs), which 
cross the entire film, are sometimes trapped in metastable states. We calculate the elastic energy 
associated with the meanderings of a TDL. The TDL behaves as an anisotropic and dispersive 
string with thermal fluctuations largely along its Burger's vector. These fluctuations also modify 
the structure factor of the vortex solid. Both effects can in principle be used to estimate the elastic 
moduli of the material. 
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PACS numbers: 74.60. Ge, 61.72.Lk, 74.72.Hs 

Recent two-sided flux decoration experiments have 
proven an effective technique to visualize and correlate 
the positions of individual flux lines on the two sides 
of BiiSr-iC aCu-iO% (BSCCO) thin superconductor films 
jj],U . This material belongs to the class of high-T c super- 
conductors (HTSC) , whose novel properties have aroused 
considerable attention in the last few years ||. Due to 
disorder and thermal fluctuations, the lattice of rigid 
lines, representing the ideal behavior of the vortices in 
clean conventional type II superconductors, is distorted. 
Flux decoration allows to quantify the wandering of the 
lines as they pass through the sample. The resulting dec- 
oration patterns also include different topological defects, 
such as grain boundaries and dislocations, which in most 
cases thread the entire film. 

Decoration experiments are typically carried out by 
cooling the sample in a small magnetic field. In this 
process, vortices rearrange themselves from a liquid-like 
state at high temperatures, to an increasingly ordered 
structure, until they freeze at a characteristic temper- 
ature Thus, the observed patterns do not repre- 
sent equilibrium configurations of lines at the low tem- 
perature where decoration is performed, but metastable 
configurations formed at this higher freezing tempera- 
ture. The ordering process upon reducing temperature 
requires the removal of various topological defects from 
the liquid state: Dislocation loops in the bulk of the sam- 
ple can shrink, while threading dislocation lines (TDLs) 
that cross the film may annihilate in pairs, or glide to the 
edges. However, the decoration images still show TDLs 
in the lattice of flux lines. The concentration of defects is 
actually quite low at the highest applied magnetic fields 
H (around 25G), but increases as H is lowered (i.e. at 
smaller vortex densities). Given the high energy cost of 
such defects, it is most likely that they are metastable 
remnants of the liquid state. (Metastable TDLs arc also 
formed during the growth of some solid films Q.) 

Generally, a good correspondence in the position of in- 
dividual vortices and topological defects is observed as 
they cross the sample. Nevertheless, differences at the 
scale of a few lattice constants occur, which indicate the 



wandering of the lines. Motivated by these observations, 
we calculate the extra energy cost associated with the 
deviations of a TDL from a straight line conformation. 
The meandering TDL behaves like an elastic string with 
a dispersive line tension which depends logarithmically 
on the wavevector of the distortion. 

b 




FIG. 1. Schematic plot of a TDL in a superconductor film. 
The Burger's vector b lies in the plane perpendicular to the 
magnetic field H. 

By comparing the experimental data with our results 
for mean square fluctuations of a TDL, it is in princi- 
ple possible to estimate the elastic moduli of the vortex 
lattice. Hence, this analysis is complementary to that of 
the hydrodynamic model of a liquid of flux lines, used so 
far to quantify these coefficients jl0| . On the other hand, 
the presence of even a single fluctuating TDL consider- 
ably modifies the density correlation functions measured 
in the decoration experiments. The contribution of the 
fluctuating TDL to the long-wavelength structure fac- 
tor is also anisotropic and involves the shear modulus, 
making it a good candidate for the determination of this 
coefficient. 

In the usual experimental set-up, the magnetic field 
H is oriented along to the z axis, perpendicular to the 
CuO-planes of the superconductor. The displacements of 
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the flux lines from a perfect triangular lattice at a point 
(r, z), are described in the continuum elastic limit, by a 
two-dimensional vector field u(r, z). The corresponding 
elastic free-energy cost is 

/d 3 r 
— [c 66 (Vm) 2 + (en - c 66 )(V • uf + c 4 4(d z u) 2 ] , 

(1) 

where V = (d x x + d y y); and en, C44, and cqq, are the 
compression, tilt, and shear elastic moduli, respectively. 
Due to the small magnetic fields involved in the experi- 
ments, non-local elasticity effects j| are expected to be 
weak, and will be neglected for simplicity. In addition, at 
the temperatures corresponding to the freezing of deco- 
ration patterns, disorder-induced effects should be small, 
and will also be ignored. 

To describe a dislocation line, it is necessary to spec- 
ify its position within the material, and to indicate its 
character (edge or screw) at each point. The latter is 
indicated by the Burger's vector b, which in the contin- 
uum limit is defined by § L du = —b, with L a closed 
circuit around the dislocation . For the TDLs in our 
problem, the Burger's vectors lie in the xy-plane, and the 
line conformations are generally described by the position 
vectors Rd{z) = R(z) + zz (see Fig. |l|). Unlike a vortex 
line, the wanderings of a TDL are highly anisotropic: In 
an infinite system with a conserved number of flux lines, 
fluctuations of the TDL are confined to the glide plane 
containing the Burger's vector and the magnetic field. 
The hopping of the TDL perpendicular to its Burger's 
vector (climb) involves the creation of vacancy and inter- 
stitial lines, as well as the potential crossing of flux lines 
@,|. These defects are very costly, making TDL climb 
unlikely, except, for instance, in the so-called supersolid 
phase, in which interstitials and vacancies are expected to 



The above expressions are obtained after integrating over 
q, with a long-wavevector cutoff A at distances of the 
order of the flux-line lattice spacing, below which the 
continuum treatment is not valid. If we also take into 
account a short- wavevector cutoff A* due to finite sam- 
ple area, the dependence of the kernels on k z has dif- 
ferent forms. For values of k z <C A*-\/min(c66, c\\)jc^ 
the logarithms in Eqs.([||4|) reduce to the constant value 
21n(A/A*), and both kernels are simply proportional to 
k 2 . In the opposite limit, if k z ^> A yj max(c66 , Ci 1 ) / C44 
all the logarithms can be approximated by the first term 
of their Taylor expansion, and A± and A11 are inde- 
pendent of k z . In between these extremes, the form of 
the kernels is globally represented through Eqs.(|^-[|). In 



proliferate P). Nevertheless, for a sample of finite extent, 
introduction and removal of flux lines from the edges may 
enable such motion. This seems to be the case in some 
of the decoration experiments where the number of flux 
lines is not the same on the two sides M. In order to be 
completely general, at this stage we allow for the possi- 
bility of transverse fluctuations in R(z), bearing in mind 
that they may be absent due to the constraints. 

We decompose the displacement field u into two parts: 
u s (r — R(z), z), which represents the singular displace- 
ment due to a TDL passing through points Rd(z) in inde- 
pendent two-dimensional planes; and u r , a regular field 
due to the couplings between the planes. By construc- 
tion, the former is the solution for a two-dimensional 
problem with the circulation constraint while the 
latter minimizes the elastic energy in Eq.rtjJ)^ and is con- 
sequently the solution to 

c 66 V 2 w r + (en - c 66 )VV ■ u r + c iA d 2 u r = ~c4id 2 z u s . 

After Fourier transforming this equation and substitut- 
ing its formal solution in Eq.(|l|), the energy cost of a 
fluctuating TDL is obtained as 7i d = W tr + AH, where 
TL str the standard energy of a straight dislocation line. 
Keeping only the lowest order terms in the small devia- 
tions R(z) of the TDL axis from the straight line, gives 
the extra energy cost of distortions as 

f dk 

AH = J^ [A±(k z )\R±(k z )\ 2 + A\\{k t )\R\\(k g )\*] . 

(2) 

Here, R(k z ) is the Fourier transform of R(z); R± and 
J2|| stand for its components perpendicular and parallel 
to the Burger's vector, respectively; and 

(3) 
(4) 



practice, the smallest wavevector k z that can be probed 
experimentally is limited by the finite thickness of the 
sample, and is ultimately constrained (by the measured 
values of en, cqq and C44) to the last two regimes. From 
Eqs.(|[^), we conclude that the TDL behaves as an elas- 
tic string with a dispersive line tension (e^ tx lnfc z ), indi- 
cating a non-local elastic energy. (A single flux line shows 
a similar dispersive behavior, as pointed out by Brandt 

Equilibrium thermal fluctuations of a TDL are cal- 
culated from Eq.(^|), assuming that one can associate 
the Boltzmann probability e -^ n / k BT ^ Q ^ s me tastable 
state. After averaging over all possible configurations of 
R(k z ), the mean square displacements are obtained as 
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(\R±(k z )\ 2 ) = 



k B TL 

A±{k z ) 



and (|i?||(A; 2 )| 2 



k B TL 



, (5) 



and given by (\R_ 



|2\l/2 



Xfim, in the same regime. 



respectively, where L is the thickness of the film. In terms 
of the function 



(\R(k z )\ 2 



8nkBTL \n-Ul 



ciiA 5 



6 2 C44fc 2 \ C44A: 2 

these quantities satisfy the simple relation 

(Mr 1 - (m 2 )- 1 = m 2 )£ (1- — 



(6) 



(7) 



Thus, even if the TDL is allowed to meander without 
constraints, its fluctuations are anisotropic, as A±(k z ) 
A (k z ) only for c\\ = c 66 . In HTSC materials, Cqq <C c u , 
so that A±(k z ) > A\\(k z ), limiting fluctuations largely to 
the glide plane. 

In real space, the width of the TDL depends on quan- 
tities such as (|i?(i)| 2 ) Cll = l/L J dk z /2n (\R(k z )\ 2 ) Cll . 
Its dependence on the film thickness L follows from 
Eq.@, as 



(\R(L)\ 2 



k B T j 1 /do 

^\ L/\2 IflnU 



if L < h, 
if L > l u 



(8) 



where we have defined l\ = by/ C44/C11, and d a is a short- 
distance cutoff along the z-axis. For length-scales below 
d a , the layered nature of the material is important. 

In order to estimate typical fluctuations for the TDL, 
we assume that cqq <C ch, and approximate Eq.(Q) by its 
leading behavior. In this limit, (\Rh | 2 ) ~ 2(|i?| 2 ) Cf56 , with 
(\R\ 2 )c 66 defined as in Eq.(||), after replacing c\\ with cqq. 
A behavior similar to Eq.(g) is obtained for (|-R| 2 )c 66 , 
with a corresponding crossover length Iq = by/ ch/cgq. 
Thus, longitudinal fluctuations of the TDL are approxi- 
mately constant for samples thinner than Z 6 , and grow 
as L/hxL for thicker samples. If allowed, transverse 
fluctuations follow from Eq.(Tj), whose leading behav- 
ior for small c 66 gives, (\R±Y) ~ (|i?| 2 ) cil . In Fig.| 
we have plotted (|i?| 2 ) Cll and 2(|i?| 2 ) Ce6 as a function 
of the thickness L. Both quantities are very sensitive to 
the elastic coefficients. We have considered the values of 
c u = 2.8 x 1CT 2 G 2 , c 44 = 8.1G 2 , and c 66 = 9.6 x 1(T 3 G 2 
reported in Ref. Q for a sample decorated in a field of 
24G, which shows a single dislocation. The Burger's vec- 
tor is equal to the average lattice spacing a a = 1/im, and 
the short-distance cutoff in the plane is taken to be of 
the same order of magnitude. The crossover lengths in- 
troduced turn out to be l\ ~ 17/im, and 1$ ~ 29^m, so 
that the experimental sample thickness (L ~ 20/im) ap- 
proximately falls into the constant regime in Fig.|^. From 
the top curve in FigJ|, we estimate (1-Ryl 2 ) 1 / 2 ~ 3/xm for 
T ~ 80K, with an uncertainty factor of about \/l0 due 
to, for instance, the uncertain values of temperature, and 
both the in-plane and perpendicular cutoffs. If uncon- 
strained, transverse fluctuations of the TDL arc smaller, 



The question mark in Fig.|2| is a reminder that once these 
fluctuations exceed a lattice spacing, proper care must 
be taken to account for constraints, and their violation 
by defects or surface effects. 
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FIG. 2. Mean-square displacements as a function of the 
thickness L, both measured in /jm, for a BSCCO sample 
decorated in a field of 24G (a = b = l^m). 

As discussed in Ref. j|, the values of cu and cqq 
measured in the experiments are about three orders of 
magnitude smaller than the theoretical predictions from 
Ginzburg-Landau theory. If we use the latter values in 
our computations, the crossover lengths become much 
shorter and TDL fluctuations are reduced by three or- 
ders of magnitude! Due to this sensitivity, analysis of 
transverse and longitudinal fluctuations of TDLs in two- 
sided decoration experiments should provide a comple- 
mentary method for determining the elastic moduli. Un- 
fortunately, the films studied so far are in the short dis- 
tance regime where details of the cutoff play a significant 
role. Experiments on thicker films are needed to probe 
the true continuum limit. 

TDLs also produce anisotropics in the flux line den- 
sity n(r,z), and the corresponding diffraction pat- 
terns. Neutron scattering studies can in principle resolve 
the full three dimensional structure factor S(q,k z ) — 
(\n(q, fc z )| 2 ), although only a few experiments are cur- 
rently available for different HTSC materials Q , ^5| . 
Two-sided decoration experiments also provide a quanti- 
tative characterization of the two-dimensional structure 
factors calculated from each surface, as well as the cor- 
relations between the two sides of the sample. 

The diffraction pattern from a vortex solid has Bragg 
peaks at the reciprocal lattice positions. Unbound dis- 
locations modify the translational correlations; a finite 
concentration of dislocation loops can drive the long 
wavelength shear modulus to zero, while maintaining the 
long-range orientational order 113]. The resulting hex- 
atic phase has diffraction rings with a 6-fold modula- 
tion, which disappears in the liquid phase. In all phases, 
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the diffuse scattering ciose to q = is dominated by 
the long wavelength density fluctuations, which are ade- 
quately described by n = V ■ it, leading to 



S(q,k z )~ (\q-u(q,k z )\ r < 



(9) 



The contribution of equilibrium density fluctuations 
(from longitudinal phonons) to Eq.(ft) has the form 110] 



S°(q,k z 



k B TLAq 2 



cn<T 



C44 



h 2 ' 

n. , 



(10) 



where A is the sample area. This contribution is clearly 
isotropic, and independent of the shear modulus in the 
solid phase. (The anisotropics of the solid and hexatic 
phases are manifested at higher orders in q.) For a sam- 
ple of finite thickness, the phonon contribution in rather 
general situations including surface and disorder effects, 
was obtained in Ref. fo|, as S°(q,L) = S° D (q)R(q, L), 
where S% D (q) is the 2D structure factor of each surface, 



while R(q, L) measures the correlations between patterns 
at the two sides of the film. 

The above results were used in Ref. (|] to determine the 
elastic moduli en and C44 of the vortex array, at different 
magnetic fields. However, the decoration images used for 
this purpose have the appearance of a solid structure with 
a finite number of topological defects. We shall demon- 
strate here that the presence of a single trapped TDL 
modifies the isotropic behavior in Eq. (|To|) . We hence- 
forth decompose the displacement field u(r) into a reg- 
ular phonon part u°, and a contribution u d = u s + u r 
from the meandering TDL described by R(z). The over- 
all elastic energy also decomposes into independent con- 
tributions H[u°, R] = H°[u°] + U d [R]. To calculate the 
average of any quantity, we integrate over smoothly vary- 
ing displacements u° , and over distinct configurations of 
the dislocation line R(z). Thus the structure factor in 
Eq.(^|) becomes a sum of phonon and dislocation parts. 

The results of Eq.(g) can be used to calculate the con- 
tribution from a fluctuating TDL, which has the form 



S d (q,k z ) = 



8nLb 2 c4 6 qi 



S(k z ) 



k B TLb 2 



"66<7||<7j_ (c 44 fcf + 2c 66 q 



1 \2 



(c n ql + C44 fe2)2 I A||(fc z ) 



A x (k z ) 



(11) 



with <7|| = q ■ b/b, and q± the component perpendicu- 
lar to the Burger's vector. The first term on the r.h.s. of 
Eq.(|il"l) corresponds to the straight TDL, and vanishes in 
the liquid state with cqq = 0. The next two terms on the 
r.h.s. result from the longitudinal and transverse fluctua- 
tions of the TDL respectively. (The latter is absent if the 
TDL is constrained to its glide plane.) The dislocation 
part is clearly anisotropic, and the anisotropy involves 
the shear modulus cqq. Thus after inverting the k z trans- 
form in Eq.(|ll|), the TDL contribution to the structure 
factors calculated from the two-sided decoration experi- 
ments can also be exploited to obtain information about 
the elastic moduli. 

In conclusion, we have calculated the energy cost of 
meanderings of a TDL in the flux lattice of a HTSC film 
such as BSCCO. Flux decoration experiments indicate 
that such metastable TDLs are indeed frequently trapped 
in thin films in the process of field cooling. We have es- 
timated the thermal fluctuations of a TDL in crossing 
the sample, as well as its contribution to the structure 
factor. Both effects can in principle be used to estimate 
the elastic moduli of the vortex solid. However, there are 
strong interactions between such defects, which need to 
be considered when a finite number of TDLs of different 
Burger's vectors are present. The generalization of the 
approach presented here to more than one TDL may pro- 
vide a better description of the experimental situation. 
From the experimental perspective, it should be possible 
to find samples with a single trapped TDL, providing a 
direct test of the theory. Other realizations of TDLs can 
be found in grown films || , and may also occur in smectic 
liquid crystals. It would be interesting to elucidate the 



similarities and distinctions between the defects in these 
systems. 
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